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D ■ There are many solutions to the Einstein field equations that demonstrate naked singu- 

Ph ' larity (NS) formation after regular evolution. It is possible, however, that such a quantum 

effect as particle creation prevents NSs from forming. We investigate the relation between 
the curvature strength and the quantum effects of NSs in a very wide class of spherical dust 
collapse. Through a perturbative calculation, we find that if the NS is very strong, the 
quantum particle creation diverges as the Cauchy horizon is approached, while if the NS is 
very weak, the creation should be finite. In the context of cosmic censorship, strong NSs will 
■») ■ be subjected to the backreaction of quantum effects and may disappear or be hidden behind 

horizons, while weak NSs will not. 
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§1. Introduction 

The cosmic censorship hypothesis (CCH) presents one of the most important 
unsolved problems in general relativity. 1 ) Its validity is often assumed in the analysis 
of physical phenomena in strong gravitational fields. There are two versions of 
this hypothesis. The weak hypothesis states that all singularities in gravitational 
collapse are hidden within black holes. This version implies the future predictability 
of the spacetime outside the event horizon. The strong hypothesis asserts that no 
singularities visible to any observer can exist. This version states that all physically 
reasonable spacetimes are globally hyperbolic. Despite several attempts, neither a 
general proof nor a precise mathematical formulation of the hypothesis has yet been 
obtained. By contrast, some solutions of the Einstein field equation with regular 
initial conditions evolving into spacetimes containing naked singularities (NSs) have 
been found. 2 )' 3 ) 

For a naked-singular spacetime to provide a counterexample of the CCH, it is 
at least necessary that the Cauchy horizon (CH) be stable. Although the CCH was 
originally stated in the classical context, CHs may be unstable due to the backreac- 
tion of quantum effects, such as particle creation, i.e., particle creation may prevent 
NSs from forming. Study of such a possibility can be traced back to the pioneering 
works of Ford and Parker 4 ) and Hiscock, Williams, and Eardley. 5 ) Ford and Parker 
considered the particle creation during the formation of a shell-crossing NS to obtain 
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a finite amount of flux. 4 ) Hiscock et al. considered the formation of a shell-focusing 
NS in the collapse of a null-dust fluid to obtain a diverging amount of flux. 5 ) Sub- 
sequently, such quantum phenomena have been studied in models of a self-similar 
dust, 6 ) -8 - 1 a self-similar null dust, 8 )' 9 ) and an analytic dust 10 ) models, for which 
the luminosities are found to diverge as negative powers of the time remaining to 
the CHs. The analytic model describes the spherical dust collapse with an analytic 
initial density profile with respect to locally Cartesian coordinates. The analyticity 
of the initial density profile and the self-similarity are incompatible in the spherically 
symmetric dust model. It is argued that the quantum radiation from a strong NS, 
such as a shell-focusing one, must diverge as the CH is approached, 6 )' 9 ) although 
this has not yet been proved. See Ref. 11) for a recent review of the quantum and 
classical radiation processes during NS formation. 

Recently, two of the present authors computed the quantum radiation in spher- 
ically symmetric self-similar collapse without specifying the collapsing matter. 8 ) It 
was found that in the generic class of self-similar spacetimes resulting in a NS, the 
luminosity of particle creation diverges as the inverse square of the time remain- 
ing to the CH. Moreover, there was another interesting result, leading us to the 
present study. In the self-similar collapse of a massless scalar field, described by the 
Roberts solution, 12 ) the luminosity remains finite at the CH. This result can be 
interpreted in terms of the curvature strength of the NS along the CH. Although 
NSs forming in generic spherically symmetric self-similar spacetimes are known 13 ) to 
satisfy the strong curvature condition (SCC) 14 ) along the CH, the NS appearing in 
the Roberts solution does not satisfy even the limiting focusing condition (LFC), 15 ) 
which is weaker than the SCC.*) The relation between the curvature strength and 
quantum effect of NSs has already been suggested in Ref. 10). It demonstrates the 
weak divergence of the quantum radiation in the so-called analytic dust model, in 
which the forming NS is known to be weak. 17 )' 18 ) However, a comprehensive un- 
derstanding of the relation between the curvature strength and the quantum effects 
of NSs has not yet been realized. The purpose here is to show how the amount 
of quantum radiation during the formation of NSs depends on such properties of 
singularities as the curvature strength. This analysis will help us obtain knowledge 
about the instability of CH, which should be predicted by a full semiclassical theory, 
taking into account the backreaction of quantum fields on gravity. In addition, it 
is shown here how the manner in which the quantized scalar fields are coupled to 
gravity changes the amount of quantum radiation. The dependence on the manner 
of coupling is important, because the CHs exhibit a semiclassical instability, caused 
by all fundamental quantum fields. 

The organization of this paper is as follows. In § 12.11 we introduce a class of 
the Lemaitre-Tolman-Bondi (LTB) solutions 19 ) that result in the formation of a NS. 
Then, the class of LTB solutions is divided into three sub-classes, depending on the 

*' Following the work of Clarke and Krolak, 16 ^ consider a geodesic (N), affinely parameterized 
by k, with tangent vector fc^, and terminating at or emanating from a singularity, where k = 0. If 
lim^^o k 2 R^k 11 k v 7^ and lim K „,o nR^k^k 1 ^ ^ 0, where R^ v is the Ricci tensor, then the SCC 
and LFC are satisfied along N, respectively. Since the quantity R^k^k" for the Roberts solution 
indeed vanishes along the CH, the NS satisfies neither the SCC nor the LFC. 
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curvature strength of the NSs in § 12.21 In § El we calculate a map of null rays 
passing near the NS. This plays crucial roles in estimating quantum emission in the 
geometrical optics approximation. In §0] the luminosity and total energy of emitted 
particles are computed. Section [3 is devoted to discussion. Throughout this paper, 
the units for which c = G = K = 1 and the signature (—,+,+,+) for spacetime 
metrics are used. 



§2. Naked singularity in spherically symmetric dust collapse 

2.1. Lemaitre-Tolman-Bondi solutions admitting a naked singularity 

The LTB solution, 19 ) which describes the collapse of a dust fluid, is written in 
a comoving coordinate system as 

T>I2 

ds 2 = -dt 2 + Trrdr 2 + R 2 (t, r)df2 2 , 

1 + fir) 

R 2 = ^ + f(r), (2-1) 



F' 



P 



8irR 2 R r 

where p is the energy density, df2 2 is the line element of a unit two-dimensional 
sphere, and the prime and dot denote the partial derivatives with respect to r and 
t, respectively. Since we are concerned with the collapse of a dust fluid, we require 
R < 0. The arbitrary functions F(r) and l+/(r) > are twice the conserved Misner- 
Sharp mass and the specific energy, respectively. In this paper, we only consider the 
case / = 0, which is called 'marginally bound collapse'. In this case, Eq. (|2-1|) is 
integrated to give 

R 3 = 9 -F(r)[t-t s (r)] 2 , (2-2) 

where t s (r) is an arbitrary function of r. The time t = t s (r) corresponds to the 
moment when a dust shell at r meets the shell-focusing singularity, which is defined 
by R = 0. It is possible to choose t s (r) = r by scaling r. LTB solutions can describe 
the formation of a shell- focusing NS from regular spacetimes. It has been shown that 
the shell-focusing singularity existing at R = with r > is completely spacelike, 20 ) 
and therefore our investigation will be confined to the singularity at (t, r) = (0,0). 

Here, we introduce a class of marginally bound LTB solutions in which the 
leading-order term of the mass function F near the regular center takes the form 

4A 3 

F{r) = 9(^ + l)3 r3M+1 + °( r3 " +1 )' ^ 

where p 0) and A (> 0) are constants.*- 1 In Appendix^! it is shown that this class 

*' This form of F is chosen so that the initial density profile can be expanded in a power 
series in R [see Eq. 12-511 1 and the singularity located at (t, r) — (0, 0) is at least locally naked (see 
Appendix[XJ. The nontrivial form of the factor and power in Eq. I2-3B are just for the simplicity of 
Eq. 12-41 . Of course, other choices of F can allow the singularity to be naked. 
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of LTB spacetimes results in the formation of a shell-focusing NS, which could be 
globally naked and therefore violate the weak version of CCH. Radial null geodesies 
are described near the center by 



^ = ±R' ~ ±\rT(t/r), (2-4) 



where 



1-^' 
3(/x+l) . 



(1-z) 



-1/3 



Here, the upper and lower signs correspond to outgoing and ingoing null geodesies, 
respectively. 

We now consider the initial regular density profile near the regular center. The 
density profile at the initial time t = t- m < is written, 



1 



1 + 2 ( 1 + ^)i + ° ((r/t '" )2) 



^to - 

Therefore, the initial density profile in terms of physical radius R oc r^ +1 / 3 takes the 
form 

Pin(R) = Po + PiR 1 + ■ ■ ■ , (2-5) 

where 

= _A_ - _L (3/i + 2)(/i + l) 3 /( 3 " +1 > 

7 - 3p + 1 ' P0 ~ 6irtl ' Pl ~ 3(3/i + l)7rA3/(3f+i) (-t in )(^+5)/(3^+i) " W 

The parameter 7 is in the region < 7 ^ 3 for p ^ 0. In particular, the analytic 
and self-similar LTB models are the cases with 7 = 2 (p = 1/6) and 7 = 3 (/i = 0), 
respectively. 

2.2. Curvature strength of the naked singularities 

The curvature strength of spacetime singularities is defined in the hope that 
weak convergence would imply the extendibility of the spacetime in a distributional 
sense. In this context, Tipler defined the strong curvature condition (SCC), 14 ) while 
Krolak defined a weaker condition, which we call the limiting focusing condition 
(LFC). 15 ) The sufficient and necessary conditions for the singularities in spherically 
symmetric spacetimes with vanishing radial pressure satisfying the LFC or SCC are 
given in simple forms. Suppose a singularity is naked, and the relation between 
the circumferential radius R and the Misner-Sharp mass m is given by 

R ~ 2y m 13 (2-7) 

near the singularity along the null geodesies terminating at or emanating from the 
NS. The constant j3 is in the region 



1/3 < 1. 



(2-8) 
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The constants yo an d are determined so that there is a positive finite limit yo = 
lim m ^o R/{2m 13 ) along those null geodesies. If the geodesic satisfies the "gravity- 
dominance condition," 18 ^ the sufficient and necessary conditions are summarized as 
the following theorem. 18 ) For the radial null geodesic terminating from or emanating 
from the NS, if and only if 1/3 < (3 < 1/2 is satisfied, neither the SCC nor the LFC 
holds, if and only if 1/2 ^ /3 < 1 is satisfied, the LFC holds but the SCC does not, 
and if and only if (3 = 1, both the SCC and the LFC hold. 

For the LTB spacetimes we consider, from Eq. (jA-2j) . we find that R oc r^ +1 
holds along the outgoing null geodesic emanating from the NS, i.e., along the CH.*) 
Also, we find that m = F(r)/2 oc r 3 ^ +1 holds along the outgoing null geodesic near 
the singularity. Therefore, the constant (3 is obtained as 

, = i±± = h±l. (2 . 9) 

3/x + 1 9 

With the above theorem, proved by Harada et a/., the NS for < 7 < 3/4 satisfies 
neither the LFC nor the SCC along the outgoing null geodesic emanating from the 
NS. The NS for 3/4 ^ 7 < 3 does not satisfy the SCC, but it does satisfy the LFC. 
The NS for 7 = 3 satisfies both the LFC and the SCC. (See also Table [I] for the 
relation between 7 and the curvature strength of NSs.) 

§3. Map of null rays passing near the naked singularity 

3.1. Local map 

The global map v = Q{u) is defined as the relation between the moments when 
a null ray leaves T~ and when it terminates at X + after passing through the regular 
center (see Fig. This global map plays crucial roles in estimating the quantum 
radiation with the geometrical optics approximation. The global map cannot be 
obtained without solving the null geodesic equation from Z~ to T + globally. 

However, the main properties of the global map should be determined by the 
behavior of null rays near the singularities, since the particle creation should be 
caused by the curvature around singularities. From this point of view, Tanaka and 
Singh proposed an alternative map, which we call the local map. 21 ^ They considered 
an observer on a comoving shell who sends ingoing null rays. These null rays are 
reflected at the regular center, and then return to the same comoving observer. A 
radial null geodesic crosses a comoving shell located at a fixed comoving radius r 
before and after the reflection at the center. Thus these null rays define a map 
between the sending time and the receiving time measured by the proper time for 
the comoving observer. (See Fig. |2]for a schematic illustration of the local map.) 

The local map is conjectured to have the same structure as the global map, 
because there are no singular features in the map between the proper time on a 



The null ray given by Eq. 1A-2II is an asymptotic solution of the null geodesic equation with 
the boundary condition t(0) = 0. In fact, Eq. 1A-2II is the earliest null ray that emanates from the 
singularity. This was proved by Christodoulou, 20 ' and that proof can be easily generalized to other 
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comoving shell at a finite distance and those measured by the null coordinates natu- 
rally defined at infinity. This conjecture has been confirmed for the self-similar dust 
model 8 - 1 ' 21 ) and the analytic dust model, 21 - 1 which is a non-self-similar spacetime. 
Two of the present authors generalized the local map to general self-similar space- 
times, and the validity of such a local map was confirmed for the self-similar Vaidya 
model. 8 ) Therefore, one can assume that the local map has the same structure as 
the global map in the LTB spacetimes. 

In this section, we calculate the local map by solving the radial null geodesic 
equation for the non-self-similar LTB spacetimes (/z > 0). The local map for the 
self-similar LTB spacetime (fi = 0) was obtained in Refs. 21) and 8). 

3.2. Outline of the calculation of local map 

To make the notation simple, we replace the coordinates r and t with r' and t' 
given by 

r' = X^r, 
t' = X^H, 

and we write r' and t' as r and t hereafter in this section. Then, the null geodesic 
equation (|2-4j) becomes 

^ = ±r^{t/r). (3-1) 
dr 

Before beginning the calculation, it is convenient to summarize the calculational 
scheme for obtaining the local map that we implement in the following subsections. 
To obtain the local map, we have to obtain solutions of Eq. (J3Tj) near r = 0, 
which correspond to radial null rays passing through the center at t = t(0) = —to 
(to > 0), and then take the limit to — * 0. We cannot expect, however, to have 
general exact solutions to this equation. Hence, here we adopt the following scheme 
to obtain the local map. First, we apply three different approximation regimes, A, 
B and C, and find three kinds of expressions, t = t A (r), t = t B (r) and t = t c (r), 
respectively. Next, we show that these three regimes have an overlapping region, 
where all three approximations are valid and the obtained approximate solutions 
can be matched with each other. Finally, we calculate the local map, that is, we 
calculate the relation between the sending time and the receiving time of the null 
ray at a comoving observer near the center. 

In § 13.31 we consider approximation regime A, in which ^ r < ^a^o is 
satisfied. Here, r?A (<S 1) is a positive constant, independent of to- I n this regime, we 
can treat the center, and then we can relate the ingoing and outgoing null rays, which 
reach the center at the same time, t = —to- In § 13.41 we consider approximation 
regime B, in which to/f/B < T < ^A^o is satisfied. Here, tjb (<C 1) is a positive 
constant. This regime can exist only when to < (ijaVb)^^^ ■ Although regime B is 
completely included in regime A, regime B enables us to have an explicit expression 
for solutions, and it is therefore essential to obtain the local map. In § 13.51 we 
consider approximation regime C, where t/r <C 1 is assumed. When we choose t/r 
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to be 0(770), where 77c (<S 1) is a sufficiently small constant, it turns out that the 
approximation is valid for to/rjc ^ f ^ Vc^- Regimes A and B obviously have an 
overlapping region. When we take the limit to — ► 0, regime B and regime C also 
have an overlapping region, and regime C is valid at finite radius. (See also Fig. 01 
for the illustration of regions, where each regime is valid.) In § 13.61 we implement 
the matching between the approximate solutions t = t B (r) and t = t c (r) in the 
overlapping region. In § 13.71 the local map is finally obtained for a comoving observer 
at a finite radius in the region for regime C. This is the generalization of the result 
of Tanaka and Singh. 21 ) 

3.3. Regime A: r < f] A tl /{1+fl) 

To find a null geodesic for ^ r < 77Aio^ 1+/ ^, we set 



to <5' 



(3-3) 



where e and 5 are constants, and C is variable for r. Here, we assume that e < 77A and 
C is 0(1). The quantity 5 is introduced for later convenience and is not necessarily 
small in regime A. From Eqs. (|3-2j) and (|3-3(l . the following relations hold: 

t = e (^+l)//^(M+l)/M 5 ( 3 . 4 ) 

r = e^ 1 )//^ 1 /^. (3-5) 

The null geodesic t = t A (r) can be expanded in e as follows: 

00 

t A (r) = -t + ^2 e (At+1)(nM+1)At ^(C), (3-6) 

n=l 

where t A (() (n = 1,2, • • • ) are functions of ( of order unity. Substituting Eq. I|3-6j) 
into Eq. ()3-l(l . we obtain the differential equations 

dt A 

^(C)=±^ +1)/ ^W-VC), (3-7) 
dt A 

-^(()=±5e- 1 T'(-5/()t A ((), (3-8) 

dt A 1 

^-(C) = ±-5^l»e- 2 T"(-5/Q(t A )\C) ± <5C^V(-<5/C)^(C), (3-9) 

and so on, where T 1 and T" denote derivatives of T with respect to its argument. 
Since we have t A (0) = — 1 , t A {0) = must be satisfied for n ^ 1. Equation (|3-7|) 
can be integrated immediately to give 

t£(C) = ±— !— + S/0 2/3 - (3-10) 

n + 1 
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With ti(() obtained above, Eq. (|3-8|) is integrated to give 



1 f< 
(C) = ^—5^ +i y^ / ^^(S/C)(l + 5/C) 2/3 dC 



n + i 



where 



h(y) 



l 



M + 1 Jo 



x 2 ^'(-l/x)(l + l/xf^dx. 



(3-11) 
(3-12) 

(3-13) 



This integration cannot be expressed in terms of elementary functions. In a similar 
way, one can write i^(C) hi the integral form 



4(0 = ± 



i 



2(M + 1) 



r 6^ +i y»tv (-6/o(i + <vo 4/3 



+( y(V-M^+i)//i^-i < ^(_ ( j/^)7 2 ^/ ( j) 



±^ t+1)(3 ^ +1)/ ^/ 3 (C/5), 



where 



l 3 (y) = r 

JO 



1 



20* + 1 



53^'(_l/a;)(l + i/ x )4/3 + x ^F>(-i/ x )i 2 ( a 



dx. (3-14) 



It should be noted that we can safely take the limit r — > in this regime, because 
we do not assume that 5 is small. 

3.4. Regime B: to/^B < t < rjAt\ 

For the approximation regime B, we additionally assume that 5 < t/b- This also 
requires that to -C 1, from Eq. (|3-4j) . The approximate solution t = t B (r) can be 
obtained by making approximation t = i A (r) with 5 -C 1. Hence, we define i^(C) as 
the function which is obtained by approximating i^(C) with (5^1. Thus, we have 



1 flOfl)//yAt+l 



(3-15) 



The approximate form of i A for can be obtained by using the asymptotic form 

of h(y) for large y in Eq. ()3T3|) given by 



h(y) = C 2 



2/i 



3(/i + 1) 2 (2// + 1) 



^+1 [l + (l/y)] 



where C2 is a constant. Except for the term C2, all terms are completely determined 
by integrating the expanded integrand. Therefore, if (C) 1S obtained as 



tf (c) = c^oh- 1 )^ 1 )/" 



2/i 



3(//+l) 2 (2/i + l 



.^+1)7^^+1 [ 1 + (,5)]. (3-16) 
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Similarly, the asymptotic form of ^3 is given by 



h{y) = c 3 + 

and hence if (C) is 



V + A* + 1 mi _ 2C 2 ' 



2/ - o7 — r-rrr 



9(/x+ l)3(2/i + 1)(3/j + 1) m 3(/x + l) 



[1 + 0(1/1/)] 



t?m = ±c^ (M+1)(3 ^ t+1)/ ^ ± 2 ^ 2 + - + 1 5 (3 M +i)/ma3m+i m , or<ni 



90 + l) 3 (2/x + l)(3// + 

+ 2C2 ^ 2 +4m+1)/m^ M + 0(5)] . (3-17) 
30 + 1) 

Then, an approximate solution is obtained from Eqs. (|3-4|) - (|3-6jl and 1)3- 15j) - (|3- 17(1 as 

i B (r) = -to + C 2 ^ +1 ±C 3 ^ +1 

1 r M+l 2 ^ r 2M+l ± V + ^ + 1 r 3M+l" 

fi + 1 30 + l) 2 (2/i + l) 90 + l) 3 (2/U + l)(3/i + l) 

x [1 + 0(5)]. (3-18) 



+ 



For later use, we have obtained explicit expression up to third order. It is straight- 
forward to compute higher orders. It should be noted that we cannot take the limit 
r -> in (I3~T%1) . 

3.5. Regime C: t/r<Cl 

Suppose i/r <C 1. Here, we introduce the condition 77c <C 1, which controls 
the order of t/r; specifically i/r is 0(r/c). In this approximation regime, we can 
expand / in t/r on the right-hand side of Eq. (|3T|) and obtain the expanded form of 
solutions. 

Let us consider the critical outgoing and ingoing null geodesies t = t±(r), which 
emanate from and terminate at the NS, i.e., t = r = 0. Then t = (r) gives the 
CH, by definition. If we assume that t/r is 0(770), these geodesies are obtained by 
expanding the null geodesic equation (|3-lj) in power series in r with the boundary 
condition t(0) = 0: 

± W H + l 30 + l) 2 (2/i + l) 

± 2 ^ / 2 + /i + / 1 -r 3 ^ +1 + 0(77 C ). (3-19) 

90 + l) 3 (2/j + l)(3// + l) v/c; v ) 

Here, the upper (lower) sign corresponds to an emanating (terminating) null ray. 
For < r < T]q^, t/r is 0(t]q) on the critical null rays, and this approximation is 

valid. That is, the approximation regime C is valid for < r < ?7 C on the critical 
null rays. 

However, we are interested in null rays that correspond to times slightly before 
the times of these critical null rays. We expand the solution as 

00 

t = t c (r) = J2^(r), (3-20) 

n=l 
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where we assume that t%(r)/r is Ofag). Substituting Eq. (jHH^ into Eq. (JHHJ and 
expanding the right-hand side, the following differential equations are obtained: 

dt c 

^- = ±r", (3-21) 
dr 

f ^3(^i)^W^9^)^ 2 ( 4 f) 2 W- < 3 ' 23 > 

These equations are integrated to yield 

t^(r) = D ± ±-^—r^\ (3-24) 
li + 1 

^^CT^- 3(, + l) 2 2 (2, + l) ^ +1 ' 

±"? (3-26) 

9(/i + l) 3 (2/i + l)(3/i + l) ' k ' 

where D± is an integration constant that comes from the integration of Eq. (13-21(1 . 
The integration constants for Eqs. (|3-22|) and (|3-23|) are set to zero. From Eqs. (|3-20|) 
and (|3-24|) - (|3-26|) . the solution takes the following form: 

t c ( r ) = D± 

± 1 r M+l 

M + l 

= - D+r^ ^ r 2 ^ +1 

^3(m+1) 3(^ + l) 2 (2/^ + 1) 

zr - Dlr^ 1 + - D+r 2 ^ ± ^ + M + 1 r 3 ' t+1 

^9(M+1) 9/x(/i + l) ± 9(^ + l) 3 (2^ + l)(3^ + l) 

+ O(^). (3-27) 

Now, we can examine in what region the solution (|3-27|) is valid. In order for 
the expansion (|3-2()|) to be valid, t„(r)/r must be 0(r/g). This is the case if 

max ^r M , — —J = 0{j]q). 

If D± is sufficiently small, the above condition implies 

D±/r] C <r< r/J/ M . 

This is the region in which approximation regime C is valid for the null geodesic 
passing through the center at t(0) = —to < 0. The relation between D± and to is 
elucidated in the following subsection. 
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3.6. Matching the approximation regimes 

Since we are interested in the null rays that are close to the critical null geodesies, 
we take the limit to — * 0. The region for approximation regime B is completely 
included in that for regime A, and the matching is trivially implemented. If we 
assume D± = O(to), the region in which regime C is valid must have an overlap 
with the region in which regime B is valid if to ^ {vaVc)^ +1 ^ ^ ■ This means that 
we can relate the integration constants D±, which appear in the regime C solution 
t = t c (r), to to, which appears in the regime B solution t = t B (r), by matching these 
two solutions in the overlapping region. 

The different expressions t = t B (r) and t = t c (r) for the null geodesic are 
obtained independently in §§ 13.41 and 13.51 We find that the solution t = t B (r) given 
by Eq. (|3- 18|) coincides with the solution t = t c (r) given by Eq. (|3-27f) at several 
lowest orders if the constant terms satisfy the following relation: 

D ± ~ -t + C 2 t^ +1 ± C z tf +1 . (3-28) 

This verifies the validity of the assumption that D± is sufficiently small and D± is 
O(to) in the limit to — > 0. This also implies that the approximation regime C is valid 
for to/i^Q <r < t]q. Here, let us see the condition on to which must be satisfied for 
the matching, to has to satisfy 

to < min ((VAVB)^^, (VAr]c)^ +1)/ ^) , (3-29) 

so that the region for the regime B and the overlapping region for the regimes B and 
C can exist. It is possible to take such a small to, because we are interested in the 
limit t -> 0. (See also Fig. H) 

3.7. Obtaining the local map 

To obtain the local map for this spacetime, we consider a comoving observer at 
r = ro, where ro satisfies ro < ??c • As time proceeds, this observer approaches the 
ingoing critical null ray t = t crit (r) and therefore enters the region to/i]c < r < 77^ [ , 
where approximation regime C is valid. Then, t±(r ) = t±(r ) are regarded as the 
sending time and receiving time of the null geodesic, where the signs ± are introduced 
to distinguish the outgoing and ingoing null rays. From Eq. (|3-27|) . the constant D± 
is specified as 

-D ± = tf\r ) - t±(r ) + O(t 77 C ), 

where t±(ro) is the moment when the observer crosses the null geodesic (|3T9|) , 
Therefore, if we take r/c to be sufficiently small, D_ (D+) can be interpreted as the 
time difference between the moments when the observer emits (receives) the null ray 
and crosses the null ray terminating at (emanating from) the NS. (See also Fig. [3) 
Since we consider a set consisting of an ingoing null ray which reaches the center 
at t = —to and an outgoing null ray which can be regarded as a reflected ray of 
the former at t = —to, we pick up both the ingoing and outgoing null rays with the 
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same value for to- When we eliminate to from Eq. (|3-28|) for both signs, we have the 
relation between D + and D- 

D_ ~D + -2C 3 (-D + f»+\ 

which can be rewritten in terms of the sending and receiving times t±(ro) at r = ro 
for sufficiently small ro as 

£(r ) ~ t c _ rit (r ) - [ffiro) - t%(r )] - 2C 3 [t c + rit (^o) - ^(ro)f +1 ■ (3-30) 

Equation (|3-3()j) is the very local map, relating the moments when the comoving 
observer locating at r = vq sends the ingoing null ray and receives the reflected 
outgoing null rays. If we revive A, the final result becomes 

iS(r„) = tl rit (r ) - [ff(r ) - i°(r„)] - 2C 3 A 3 [tf\r ) - t%{r f^ 1 . (3-31) 

Now, note that the comoving observer must be in the region ioA?c % r o ^5 ffc^A -1 /^. 
This implies that the asymptotic structure of the local map in the limit to — ¥ 0, 
therefore the main feature of the global map, is determined by only the behavior of 
the null rays in the small but finite region < r < n^/^ X^ 1 ^ 1 . 

§4. Luminosity and energy of particle creation 

We consider the massless scalar fields coupled to the scalar curvature R as 

(0-^)0 = 0, (4-1) 

where £ is an arbitrary real constant. In particular, the scalar fields with £ = and 
£ = 1/6 are minimally and conformally coupled ones, respectively. The stress-energy 
tensor of the scalar field in an asymptotically flat region is given by 

where is a flat metric. To calculate the luminosity of emitted particles we need 
an estimate of the vacuum expectation value of the above stress-energy tensor. A 
suitable regularization is required in the calculation of the vacuum expectation value, 
because the stress-energy tensor is quadratic in fields at the same point. The regu- 
larization for minimally and conformally coupled scalar fields is obtained in Ref. 4) 
using the point- splitting regularization scheme. This can easily be generalized to the 
above more general scalar fields to give 



H(f) 2 +H)f 



(4-2) 



where the prime denotes differentiation with respect to u, and I (= 0, 1, ■ • • ) and m (= 
0, ±1, • • • , ±Z) are angular momentum and magnetic quantum numbers, respectively. 
The luminosity is the sum of all these modes, but it diverges, because the luminosity 
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given in (j4-2j) is independent of (l,m). Such a divergence is due to the fact that we 
ignored the back scattering by the potential barrier, which will reduce the emission 
for highly rotational modes. Hereafter, we omit the quantum numbers (l,m), and 
it should simply be kept in mind that the above expression holds only for small I. 
The total energy of emitted particles is estimated by integrating the luminosity with 
respect to u: 

E®(u) = [ U L^(u')du'. 



4.1. Quantum radiation in the non- self- similar LTB spacetimes: < 7 < 3 

As described in § EJ we assume that the local map and the global map have the 
same structure. This means that from Eq. ()3-31j) . the asymptotic form of the global 
map will take the form 

Q(u) ~ vq - A{uq - u) - AgX 3 (u - u) 3fl+1 , 

where u = uq and v = vq are the CH and the ingoing null ray that terminate at the 
NS, respectively, and A and g are constants. Keeping in mind that \x = (3 — r y)/{S r y), 
one can calculate the luminosity as 



+ 



r 7 7 ~ 15 
36( 7 - 2) 



27( 7 -2)(3- 7 ) 2(3- 7 )/ 7( _ n) -4 +6/7 

47T7 4 



+ o(( Uo -n)- 5+9 ^), (4-3) 

where we have defined the "frequency" of the NSs as 

Ws = A 1 ^ = A 37/(3 ^ 7) . (4-4) 

In Appendix [0 we show that the frequency (|4-4f) is identical to that defined 
in Ref. 10). Depending on whether or not £ = 1/6, the leading-order term of the 
luminosity changes. Let us examine the time dependence of the luminosity and the 
total energy of emitted particles in detail for the cases £ 7^ 1/6 and £ = 1/6. 

In the case £ 7^ 1/6, the first term in Eq. (|4-3|) dominates, except in the special 
case 7 = 3/2. For < 7 ^ 1, the leading term vanishes as u — > uq. For 1 < 7 < 3 and 
7 7^ 3/2, the luminosity diverges as the CH is approached as a negative power of the 
time remaining to the CH. The special case 7 = 3/2, in which the factor of the first 
term in Eq. (|4-3|) vanishes, is divided into two cases, depending on whether £ 7^ 1/4 
or £ = 1/4. If £ 7^ 1/4, the second term in Eq. (|4-3|) survives as a finite constant. If 
£ = 1/4, the factor of the second term also vanishes. Therefore, higher-order terms 
contribute to the luminosity, which is finite for 7 = 3/2. Therefore, in both cases 
that 7 = 3/2, the luminosity remains finite at the CH. Now, let us examine the total 
energy of emitted particles. For 3/2 < 7 < 3, the leading term is 

^^rf-^>o " u)-^'\ (4-5) 
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and therefore the total energy diverges as the CH is approached. In particular, in the 
case of an analytic LTB solution (7 = 2), the energy diverges as (uq — u)~ 1 / 2 , which 
coincides with the result in Ref. 10). For < 7 < 3/2, the total energy remains 
finite as the CH is approached. In the special case 7 = 3/2, the energy also remains 
finite. [See also Table Ufa).] 

In the case of the conformally coupled scalar field, which is defined by £ = 1/6, 
the second term in Eq. (|4-3j) dominates. It is found that in the case < 7 ^ 3/2, the 
luminosity remains finite at most. While, for 3/2 < 7 < 3, the luminosity diverges, 
due to the negative power of the time remaining to the CH. Let us examine the total 
energy of the emitted particles. For 2 < 7 < 3, we have 

which diverges as the CH is approached. For 7 = 2, we have 

2507T 

and therefore the energy diverges logarithmically. This too is identical to the result 
in Ref. 10). If < 7 < 2, the energy remains finite, at most. [See also Table Hfb).] 

It follows that the radiation resulting from a conformally coupled scalar field 
is milder than that resulting from non-conformally coupled scalar fields for a given 
value of 7. Such a consequence appears to result from the fact that the coupling of the 
conformal scalar field to gravity is weaker than that of the other scalar fields. In this 
relation, it should be pointed out that a conformally coupled scalar field must have 
non-zero and finite mass to be created in the early universe, while non-conformally 
coupled scalar particles are created regardless of their mass. 22 ) 

4.2. Quantum radiation in a self-similar LTB spacetime: 7 = 3 

The global map for self-similar LTB spacetimes ending in NS formation is calcu- 
lated analytically in Ref. 6). Subsequently, its main property was re-produced with 
the local-map method. 8 )' 21 ) The global map for null rays passing near the CH is 
given by 

G{u) ~ v - B(u - u) a [1 + q(u -u) + ((n - u) 2 )] , (4-8) 

where a, B, and q are constants. The terms in the square brackets in Eq. (J4-8|) 
together form an analytic function of (uq — u),*^ The constant a depends only on 
the parameter A in Eq. I|2-3|) . and it has been shown to be greater than unity for the 



*' In previous works on particle creation during NS formation in the self-similar LTB spacetime, 
only the constant term in the square brackets in Eq. 14-81 was considered. 6 '' 8) ' 21) This was sufficient 
to obtain results. It is easy to calculate the higher-order terms using the local-map method and 
to show that they constitute an analytic function near the CH. It is possible, however, that the 
emergence of the scale q in Eq. 14-80 from such a scale invariant spacetime as self-similar LTB 
solution indicates the breakdown of the local-map method. It is not clear to what extent the local- 
map method is valid. 
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region of A in which the singularity is naked. 6 - 1 Using the global map ([4-8)1 . 
compute the luminosity and energy of the particle creation as 



we can 




(4-10) 



(4-9) 



The first term in both Eqs. 1)4-9)1 and 1)4-10)1 dominates, except in the special case 
that £ = (a + 1)/12. Therefore, the luminosity and energy generically diverge as 
the inverse square and the inverse of the time remaining to the CH, respectively. 
However, in the special case £ = (a + 1)/12, each second term in Eqs. 1)4-9)1 and 
l)4TU)l dominates so that the power and energy diverge inversely and logarithmically, 
respectively. This case, however, should be regarded as an accidental cancellation 
in the sense that the coupling constant £ happens to be a special value a, which is 
determined by the details of the collapse. [See also Tables [IJa) and (b).] 



In this paper we have considered particle creation during the formation of shell- 
focusing NSs in a wide class of spherical dust collapse, which is described by the 
marginally bound LTB solutions. Each solution has a different initial density profile, 
and the resulting NSs have a variety of curvature strengths along the CHs. The 
luminosity and energy of particle creation have been estimated for each LTB solution 
and each scalar field that couples to scalar curvature in the linear form. The results 
are summarized in Tables QJa) and (b). 

We first mention the validity of the approximations which have been made in 
this article. Next, we discuss the relations between the quantum radiation and the 
curvature strength of the NSs and also the manner in which the scalar fields couple 
to gravity. Finally, we discuss the implications of the present result to the CCH. 

The analysis has been based on three assumptions: the validity of the local-map 
method, the geometrical optics approximation, and quantum field theory in curved 
spacetime. The validity of each approximation seems to be open to debate. (See 
Ref. 10) for discussion of the geometric optics approximation and Refs. 28) and 8) 
for quantum field theory in curved spacetime.) Here, we focus on the local-map 
method, on which our entire analysis is based. The point is that the crucial factor 
of particle creation, the redshift of particles, must be determined by the geometry 
near the NS, while in Hawking radiation, the redshift is determined by the event 
horizon, which is a global object. It is unlikely that the global map has a structure 
that differs from that of the local map, since there is no singular feature in the map 
between the moments on the comoving observer at a finite distance and that of the 
null coordinates naturally defined at infinity. Indeed, in the models of the self-similar 
LTB, 8 - 1 ' 21 - 1 the analytic LTB 21 ) and the self-similar Vaidya, 8 ) the local-map method 
yields the correct results, which are obtained with the global map. For this reason, 
we have assumed the validity of the local-map method. 

From the results, it is found that following statements hold for the generic naked- 
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singular LTB spacetimes: the SCC along the CH is a sufficient condition for the lu- 
minosity and energy of the created scalar particles to diverge as the CH is approached; 
while, violation of the LFC is a sufficient condition for the luminosity and energy to 
be finite at the CH; if the NS does not satisfy the SCC but does satisfy the LFC, the 
luminosity and energy can be either divergent or finite. We only consider dust col- 
lapse for simplicity. However, the above statements regarding the curvature strength 
and the quantum radiation should be independent of the collapsing matter, because 
particle creation is a purely kinematic phenomenon and not directly related to the 
Einstein field equations. Therefore, we conjecture that the above statements hold 
for spherically symmetric collapsing spacetimes with any kind of collapsing matter. 
Of course, its validity should be verified or examined with known solutions ending 
in NS formation. The self-similar models which have been examined to this time, 
the collapse of a null-dust fluid, 8 )' 9 )' 23 ) a massless scalar field, 8 )' 12 ) and a perfect 
fluid, 8 )' 24 ) support this conjecture. There are many examples to be investigated: 
NS formation in counter-rotating particles, 18 )' 25 ) non-self-similar null dust, 23 ) null 
strange quark matter, 26 ) various matter fields in the critical collapse 27 ) and so on. 
Here, we also mention the manner in which scalar fields couple to gravity. Although 
the quantum radiation due to the conformally coupled scalar field is less than that 
of other scalar fields, including the minimally coupled scalar field, the dependence 
of the amount of radiation on the manner of coupling is not so drastic as to modify 
the above statements. 

Next, we consider the implications of our results for the CCH. The diverging 
radiation from strong NSs corresponds to an instability of the strong NS formation. 
The system will enter into a phase where the backreaction from the quantum field 
to spacetime plays an important role. While, the finite radiation from the weak NSs 
corresponds to a stability of the weak NS formation. This is striking, because the 
weak NSs seem to need another mechanism if they are to hide behind horizons. Of 
course, we cannot dismiss the possibility that the effect of backreaction suppresses 
the quantum radiation and that strong NSs appear, all things considered. 

In the present analysis, we did not find a necessary and sufficient condition on 
the curvature strength of NSs for the quantum radiation to be divergent or finite. 
We believe that a new definition of the strength of (naked) singularities should be 
proposed from the viewpoint of the behavior of quantum fields on spacetimes rather 
than the viewpoint of the behavior of classical particles. Such a philosophy can also 
be seen in the wave-probe approach to NSs, 29 ) which is based on the theory of 
dynamics in non-globally hyperbolic spacetimes developed first by Wald. 30 ) With 
regard to this point, there is room for further investigation. 
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Appendix A 

Nakedness of the Singularity 

In order to determine whether or not the singularity is naked, we investigate the 
future-directed outgoing null geodesies emanating from the singularity at (t, r) = 
(0, 0). We find the asymptotic solutions that obey a power law near the center 31 ) in 
the form 

t ~ X rP, (A-l) 

where Xq > and p ^ 1 are constants. The latter condition is due to the fact that 
the orbit of the shell- focusing singularity is t = t s (r) = r. After some straightforward 
calculations, for \x > we find the asymptotic solution 

^— r^ +1 . (A-2) 
H + 1 

With Eq. (|A-2(I and the fact that the apparent horizon, which is defined by F = R, 
behaves as t = t a h(r) = r — 2F(r)/3 — r for \x > near the center, the singularity 
is at least locally naked. In the case of the self-similar case {jx = 0), a similar 
treatment is possible, and the singularity is known to be naked for small values of 
A. 31 ) We consider the situation in which the collapsing dust ball is attached to 
an outer vacuum region at a comoving radius r = [constant], within which the null 
ray (|A-2|) is outside the apparent horizon. Then, the singularity is globally naked, 
and the weak version of CCH is violated. 

Appendix B 

Frequency of the Singularity 

It is helpful to compare the gauge used in this paper with one used in much 
of the literature. This comparison shows that lo s defined in § 0] coincides with the 
characteristic frequency of the singularity introduced in Ref. 10), up to a numerical 
factor. Let us denote the comoving coordinates by (t, f), in which f is chosen to 
coincide with the physical radius R at the initial regular epoch t = 0, i.e., we have 
i?(0, f) = r. We assume that the mass function F(f) can be expanded near the 
regular center as 

F(f) = F x r a + F 2 f b + • • • , 

where a and b are constants satisfying a < b. Then the initial density profile can be 
written 



P (o,f) = ^- 3 + ^- 3 + 



(B-l) 
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Comparing Eq. (jB-l|) with Eqs. (|2-5j) and (|2-6j) . we obtain the powers and coefficients 
of F{f) as 

o h 3(3^ + 2) 
a = 6, o ~ 



3/i + l 1 

8 ( M + 1)3/(3/^+1) 



9t? n ' * 9A 3 /( 3 M+l)(-t in )(9M+5)/(3/,+l)- 

It is found that the power 6 satisfies 3 < b ^ 6 for // ^ 0. In Ref. 10), Harada ef al. 
determined the characteristic frequency of the naked singularity in the analytic model 
(fi = 1/6) using physical consideration. It is easy to repeat their computations for 
general values of /i > 0. One possible quantity, composed only of F\ and which is 
independent of the choice of the initial time slice and has the dimension of frequency 
is 

F (9M+5)/(6m) (_ F2 )-(3 M +1)/(3m) 

(9 M +5)/(3/ t ) / q \ (3 M +1)/(3 M ) 



/ 9 \ ty/i+bJA-w / q \ (isn+xj/w 



where we have used Eq. (|4-4|l . In terms of 7, the quantity ()B-2|) can be written as 
follows: 

f2 7 (\) = F (^yw-^\-F 2 )- 3 /^\ 

In the case of the analytic LTB model (7 = 2), we have 

^ 2 (A) = F 1 13/2 (-F 2 )- 3 , 



which coincides with the frequency defined in Ref. 10), up to a numerical factor. 
This shows that it is valid to define uj s = X 1 ^ 1 as the frequency of a singularity. In 
the self-similar LTB solution (/z = 0), such a quantity does not exist, because of the 
scale-invariant nature of self-similar spacetimes. 
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Table I. The relations among the curvature strength of naked singularities, the luminosity, and 
the energy of scalar fields near the Cauchy horizons (a) for the non-conformally coupled scalar 
fields (£ =^ 1/6) and (b) for the conformally coupled scalar field (£ = 1/6). The constant 7 
parameterizes the initial density profile of the dust fluid. SCC implies LFC. 
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, . r = U 
Singularity t > I 




i 



Fig. 1. A possible causal structure of the LTB spacetimes considered in this article. A null singu- 
larity occurs at the spacetime point (uo,Vo), and it is visible from I + , where (u, v) is a suitable 
double null coordinate system. An outgoing null ray u = [constant] can be traced backward in 
time from T + to I - . It turns out to be an ingoing null ray v — Q{u) after passing the regular 
center, located at r = 0, with t < 0. The outgoing null ray u — uo and the ingoing null ray 
v — vq represent the CH and the null ray that terminates at the NS, respectively. 
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Fig. 2. A schematic spacetime diagram of the region filled with dust, with the illustration of the 
local map defined in § [21 A pair consisting of an ingoing null ray t_(r) and an outgoing null 
ray t+(r) is depicted (solid curve). It passes near the NS, located at (t,r) — (0,0) . Null 
rays terminating at and emanating from the NS, t_ (r) and t+(r), are also depicted (dashed 
curves). The latter is the CH. A comoving observer is located at r — ro — [constant]. The local 
map is defined as the relation between t-(ro) and t+(ro). 
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> r 



ri A tl^ +1) 
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C 



-> r 



1/H 

Vc 



Fig. 3. A schematic illustration of the regions in which the regimes A, B and C are valid. The 



1/(1+/*) 



to/m < r < rj A t 



regions corresponding to the regimes A, B and C are ^ r < r]At 
and D±/r]c ~ to/f]c J$ r < Vc^' respectively, where 7?x <C 1 (X = A,B,C) is a con- 
stant independent of to. The region for B is included in that for A. The regions for B 
and C exist if to < (^a^b)' m+1 ^ m and to < rjQ are satisfied, respectively. The re- 

gions for B and C overlap if to < (?7a'7c) < ' m+1 '' /m - All these conditions are satisfied for 
to < min[(77A??B) (M+1)/M , (r?AJ7c) (M+1)/M ]- In the limit to -> 0, the regions for A and B shrink to 
zero, but the region for C, where the comoving observer should be located, remains finite. 



